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The  three-layer  balanced  axisymmetric  tropical  cyclone 
model  developed  by  Ooyama  < 1969)  is  generalized  by  relaxing 
the  balance  approximation  and  replacing  the  incompressible 
fluid  layers  with  compressible  fluid  layers  in  isentropic 
coordinates.  In  numerical  simulations  the  effects  of 

latitude  are  non-linear;  storms  that  intensify  the  fastest 

u  ft  i 

< at  i0y,  N>  do  not  intensify  the  most.  When  latitude  is 
changed  with  other  parameters  held  constant,  storms  between 

■  l ( 

20  and  30^N  intensify  the  most.  When  the  latitudinal 

change  in  sea  surface  temperature  is  considered,  storms 

^  ‘,d 

between  15  and  25/N  intensify  the  most.  The  latter  agrees 
with  climatological  trends  observed  in  the  western  Pacific 
during  1980-1982.  The  effect  of  latitude  is  the  dominant 

control  of  minimum  surface  pressure  attained  for  lower 

,  >  ) 

latitudes  <10^N);  whereas  northward,  the  effects  of 
latitude  and  sea  surface  temperature  are  dominant.  The 
size  of  storms  is  largely  dependent  upon  latitude,  where 
larger  storms  occur  at  more  northerly  latitudes. 

Latitude  influences  the  development  of  tropical 
cyclone  growth  by  controlling  the  magnitude  and  radius  of 


?h 

i 
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maximum  updraft.  At  lower  latitudes#  the  zone  of  maximum 
updraft  occurs  closer  to  the  center  of  the  storm  producing 
1  aerially  smaller  storms#  and  the  rate  of  heating  is  greater 

than  for  more  northerly  storms.  Both  factors#  combined 
with  the  increased  heating  efficiency  caused  by  higher 
|  inertial  stability  values#  allow  the  storm  to  intensify 

rapidly. ^However#  because  the  rapid#  more  efficient  heating 
quickly  decreases  the  vertical  instability#  the  period  of 
intensif ipat ion  is  shorter  than  at  higher  latitudes# 
producing!  only  moderate  to  weak  intensity  storms. 
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Abstract 


PICKLE#  JOHN  D.  Tha  Effect  of  Latitude  on  the  Development 
of  Tropical  Cyclones  (Under  the  direction  of  Mark  DeMaria>. 

The  three-layer  balanced  axisymmetr ic  tropical  cyclone 
model  developed  by  Ooyama  (1969)  is  generalised  by  relaxing 
the  balance  approximation  and  replacing  the  incompressible 
fluid  layers  with  compressible  fluid  layers  in  isentropic 
coordinates.  The  development  of  a  cumulus  parameter isat ion 
scheme  is  simplified  in  isentropic  coordinates.  The 
primitive  equations  are  solved  using  the  spectral 
(Galerkin)  method  with  Bessel  function  basis  functions.  In 
numerical  situations#  an  increase  in  sea  surface 
temperature  or  mid-level  moisture  Increased  the  amount  of 
moisture  available  for  condensation  per  unit  volume  of 
ascending  air.  This  resulted  in  a  near-linear  decrease  in 
the  minimum  surface  pressure  attained#  exponential 
increases  in  the  maximum  rate  of  intensification#  and 
approximately  no  change  in  the  size  of  storms  (based  on  the 
radius  of  30-knot  winds)  at  maximum  intensity. 

The  effects  of  latitude  are  non-linear;  storms  that 
intensify  the  fastest  (at  10°N)  do  not  intensify  the  most. 
Nhen  latitude  is  changed  with  other  parameters  held 
constant#  storms  between  20°  and  30°N  intensify  the  most. 
When  the  latitudinal  change  in  sea  surface  temperature  is 
considered#  storms  between  15°  and  25°N  intensify  the  most. 
The  latter  agrees  with  climatological  trends  observed  in 


the  waitirn  Pacific  during  1980-1982.  Tha  affact  of 
latituda  is  tha  dominant  control  of  minimum  surface 


pressure  attained  for  lower  latitudes  (10°N>;  whereas 
northward*  the  effects  of  latitude  and  sea  surface 
temperature  are  dominant.  The  size  of  storms  is  largely 
dependent  upon  latitude*  where  larger  storms  occur  at  more 
northerly  latitudes. 

Latitude  influences  the  development  of  tropical 
cyclone  growth  by  controlling  the  magnitude  and  radius  of 
maximum  updraft.  At  lower  latitudes*  the  zone  of  maximum 
updraft  occurs  closer  to  the  center  of  the  storm  producing 
aerially  smaller  storms  and  the  rate  of  heating  is  greater 
than  for  aore  northerly  storms.  Both  factors*  combined 
with  the  increased  heating  efficiency  caused  by  higher 
inertial  stability  values*  allow  the  storm  to  intensify 
rapidly.  However*  because  the  rapid*  more  efficient  heating 
quickly  decreases  the  vertical  instability*  the  period  of 
intensification  is  shorter  than  at  higher  latitudes* 
producing  only  moderate  to  weak  intensity  storms. 
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Tropical  cyclonas  with  thair  daatruct ivonaaa  and 
complexity  have  been  the  focue  of  numeroue  theoretical  and 
observat ional  etudiee.  The  deetruct iveneee  ia  a  result  of 
the  intense  low-level  winds  and  the  torrential  rains.  The 
complexity  arises  from  interactions  between  scales  of 
motion  that  are  dynamically  quite  different#  i.e.  synoptic 
scale  and  cumulus  convection  scale.  Ooyama  (1964)  and 
Charney  and  Eliassen  <1964)  have  shown  that  circulations 
of  these  distinct  scales  could  enhance  each  other  by 
cooperative  interaction#  often  referred  to  as  Conditional 
Instability  of  the  Second  Kind  <CISK>. 

Merrill  ( 1985)  summarised  the  evolving  conceptual 
understanding  of  tropical  cyclone  structure  from  1900  to 
1980.  From  1900  to  1957  the  large-scale  structure  of  the 
storms  were  qualitatively  studied  from  ground-based 
stations.  The  advent  of  flying  instrumented  aircraft  into 
tropical  cyclones  and  the  use  of  computers  from  1958  to 
1963  brought  about  a  concentrated  study  of  the  core  of  the 
hurricane.  Most  of  the  observational  data  was  within  100 
km  of  the  storm  center  as  a  result  of  the  limited  flying 
range  of  the  aircraft.  Between  1964  and  1975#  numerical 
modeling  of  tropical  cyclones  rapidly  developed  as  a  result 
of  the  concept  of  CISK  and  Project  8T0RMFURY#  a  plan  to 
reduce  the  maximum  winds  in  a  hurricane.  During  this  time 
research  on  environmental  factors  essentially  stopped  and 


the  study  of  large-scale  aspects  of  hurricanes  was  outside 


the  mainstream  of  hurricane  research.  Interest  in  the 
large-scale  aspects  of  tropical  cyclones  renewed  as  a 
result  of  the  developments  of  more  advanced  weather 
reconnaissance  airplanes*  operational  geosynchronous 
satellites  with  infrared  sensors*  and  the  combining  of 
observations  from  various  hurricanes  into  oomposites 
during  1976  and  1980.  One  of  the  major  limitations  of  the 
previous  research  is  that  the  primary  emphasis  has  been  on 
the  detailed  structure  near  the  storm  center  rather  than 
the  full  horizontal  extent  of  the  storm  circulation. 

Several  studies  have  dealt  with  quantifying  the  size 
of  tropical  cyclones  and  analyzing  environmental  factors 
that  may  influence  the  storm  size  (e.g. *  Merrill  1982; 
Weatherford*  1985;  Brand*  1972;  and  McBride  and  Zehr* 
1981) *  and  a  few  papers  have  noted  the  effect  of  sea 
surface  temperature  or  latitude  on  the  size  of  the 
hurricane  in  numerical  models  < Ooyama*  1969;  Yamasaki* 
1968).  The  purpose  of  this  research  is  to  develop  a 
numerical  model  with  maximum  simplification  that  can  be 
used  to  analyze  the  development  of  tropical  cyclones.  The 
model  will  then  be  used  to  investigate  factors  which  affect 
the  storm  size.  The  results  will  be  compared  to 
climatology  data  of  Air  Force  weather  reconnaissance 
flights  through  tropical  cyclones  in  the  western  Pacific 
during  1980  to  1982  (Weatherford*  1985). , 
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An  axisymmetr ic  numerical  modal  was  chosen  because 
many  of  the  physical  aspects  of  tropical  cyclones  are  not 
critically  dependent  on  azimuthal  variations  (Anthes#  1982) 
and  the  computational  memory  and  run  time  of  an 
axisymmetric  model  are  significantly  less  than  for  a 
three-dimensional  numerical  model.  One  of  the  first 
axisymmetric  numerical  models  that  produced  a  circulation 
and  evolution  that  resembled  an  average  tropical  cyclone 
mas  developed  by  Ooyama  (1969).  The  model  to  be  used  in 
this  study  is  similar  to  Ooyama' s  model#  but  is  slightly 
more  general.  Both  models  contain  a  relatively  simple 
cumulus  parameterization#  the  minimum  vertical  resolution 
(three  layers)#  and  the  simplest  geometry  (f-plane)  which 
can  be  used  for  tropical  cyclone  simulations.  The  two 
primary  modifications  are  that  compressible  fluid  layers 
with  isentropic  coordinates  rather  than  incompressible 
fluid  layers  are  used  and  that  primitive  equations  rather 
than  the  balance  equations  are  used. 

Another  modification  to  Ooyama' s  model  is  that  the 
governing  equations  are  solved  by  the  spectral  method 


rather 

than 

the 

finite  difference 

method. 

8pectral 

methods 

are 

much 

more  accurate 

per 

degree  of 

freedom 

compared 

to 

finite  difference 

methods#  and 

spectral 

methods  can  reduce  computational  dispersion  caused  by  an 


artificial  phase  speed  for  waves  with  wavelengths  close  to 
twice  the  grid  spacing  in  a  finite  difference  model 


< Machenhauer*  1979;  DeMaria*  1983).  In  meteorology*  the 
spectral  method  has  been  used  in  global  atmospheric 
modeling  since  the  early  1970s*  but  the  method  has  been 
employed  in  modeling  tropical  cyclones  since  1983 
(DeMaria*  1983). 

In  the  spectral  method  the  spatial  dependence  of  the 
dependent  variables  is  represented  by  a  finite  series  of 
basis  functions  which  are  orthogonal  with  respect  to  some 
inner  product  and  satisfy  the  same  boundary  conditions  as 
the  dependent  variables.  The  equations  for  the  time 
dependent  amplitudes  are  found  by  substituting  the  series 
expansions  into  the  governing  equations  and  then  taking 
the  inner  product  of  the  equations  with  each  of  the  basis 
functions.  The  governing  equations  are  then  transformed 
into  a  series  of  ordinary  differential  equations  for  the 
expansion  coefficients  (often  referred  to  as  spectra). 
When  these  equations  are  solved  using  time  differencing  it 
is  necessary  to  evaluate  the  transform  of  the  nonlinear 
terms  at  each  time  step.  Using  the  transform  method 
< Orszag*  1970)  the  nonlinear  terms  are  computed  by  first 
transforming  the  dependent  variables  to  physical  space  on 
some  specific  grid.  The  nonlinear  products  are  then 
calculated  at  the  grid  points  and  the  inverse  transform  is 
evaluated  using  an  appropriate  numerical  quadrature  rule. 
The  transform  method  makes  the  inclusion  of  physical 
processes  straightforward  since  the  dependent  variables 


are  calculated  in  physical  space  at  each  time  step 


Machenhauer  <  1979>#  DeMaria  (  1983)#  and  Schubert  and 
DeMaria  ( 1984)  provide  a  more  detailed  history  and 
description  of  the  spectral  method. 
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The  governing  equations  for  a  rotating*  axiayausetric 
hydrostatic  fluid  in  iaentropio*  cylindrical  coordinates 
can  be  written  as: 


du 

du  .  du 

v2 

dM 

—  + 

u —  +  e — 

—  —  — 

fv  +  — 

=  p 

C  1> 

at 

dr  d0 

r 

dr 

u 

dv 

dv  ,3v 

uv 

—  + 

u —  +  0 — 

+  —  + 

fu 

«  F 

(2) 

at 

dr  d0 

r 

V 

dM 

~  m 

n 

<3> 

ae 

a* 

1  d(r<ru> 

d<0<r> 

—  + 

- +  - 

«  0 

(4) 

at 

r  dr 

d0 

d0 

Q 

—  m 

0  ■  - 

<5> 

dt 

n 

o  *  — 


n  «  cp<p/p0> 


N 

where: 


o  t  +  <t>  *  en  +  $ 

p 


u  *  radial  component  of  the  wind 
v  ■  tangential  component  of  the  wind 
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6  “  potential  temperature 


•  ae 

9  ■  —  ■  vertical  velocity 
dt 

f  *  Coriolis  parameter 
t  «  time 

r  *  radial  distance  from  the  storm  center 

&  ®*v  “  friction  terms 

o  -  pseudo-density 

M  »  Montgomery  stream  function 

n  =  Exner  function 

p  =  pressure  <  p  *  lOOOmb) 

cp  =  specific  heat  of  dry  air 

K  =  R .  /  c 
d  P 

T  *  absolute  temperature 

<t>  a  geopotential  *  g*  (9s  acceleration  of 

<  gravity 

(  s  *  vertical  height 

Q  a  heating  term 
dQ 

Q  =  —  =  rate  of  heating 
dt 

Assuming  Q  and  F  are  known  or  parameterized*  then 
equations  (l)-<6)  are  six  equations  with  independent 
variables  r*  6*  and  t;  and  dependent  variables  u*  v*  M*  o  , 
n,  and  p. 


7 


I 


2.2  Vertical  Discret irat ion 

For  upper  and  lowar  boundary  conditions*  aasuma  that 
p  *  Pjj  *  constant  1 


0  =  0^  =  constant 


upper  boundary  <  9  > 


0=0^  =  constant  j  lower  boundary  (10) 


To  vertically  discretise  equations  <l>-<8>*  consider  a 
fluid  which  has  three  layers  of  constant  potential 
temperature  as  shown  in  Figure  1. 


At  r  =  0.  M  =  0 II  so 


M0  •  9on. 


<  11) 


Integrating  equation  <3>  across  the  interfaces  between 
layers  0-1  and  1-2  gives: 


**1  -  »o  +  <*i 

’  Vni/2 

(  12) 

M2  =  +  <02 

“  ffl>n3/2 

(  13) 

The  momentum  equations* 

< 1) — <  2)  are 

more  easily 

discretised  when  in  flux  form.  Therefore*  multiply  < 1>  by 
a  and  (4)  by  u  and  adding  gives: 


d<ua>  1  d(u^r<r>  d<  u<t0  >  <rv2  dM 

-  + - + - fff  v  +  —  =  a  F„  (  14  > 

dt  r  dr  d9  r  dr 


u 


8 


Layer 

2 

e2' 

«2' 

u2' 

v2 

Layer 

1 

er 

ur 

V1 

Layer 

0 

e0' 

M0' 

U0' 

v0 

n5/2'  P5/2'  WT  = 
n3/2'  P3/2'  w3/2 
n  1/2'  Pl/2'  Wl/2 

nB*  pb*  «.  *  o 


o 


Figure  1  Schematic  representation  of  the  three  layers  of 
the  numerical  simulation  and  the  notation  of  the 
basic  dependent  variables. 
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and  multiplying  <2>  by  a  and  (4)  by  v  and  adding  gives: 


d<  v<7>  1  d(uvra)  d<  v<r0)  <ruv 

dt  r  dr  d0  r 


<  15) 


Since  9  is  constant  in  each  layer*  then*  from  <3>*  M 


will  also  be  vertically  constant  in  each  layer.  Since  the 


pressure  gradient  force  is  constant  in  the  vertical  for 


each  layer*  it  is  also  assumed  that  u  and  v  are  also 


vertically  constant  in  each  layer. 


Vertically  integrating  < 14)  and  (15)  across  layers  0* 


1*  and  2  gives: 


>0  Ap0)  1  a<u02rAp0>  +  -  AP0V05 

- - + - +  w,  ._u0  -  w,  , - 


r  dr 


1/20  1/21 


(  16a) 


"  Ap0fv0  +  Ap0“"  *  Ap0Fu0 


d(  u^Ap^)  1  d<u12rAp1>  +  -  + 

—  "*■  W-3  J  ->U0  —  w  ,  .  _U, 


r  dr 


3/21  3/22  1/20 


<  16b) 


Ap  i  v  ] 


on. 

+  W1/2U1 - - -  "  APl**!  +  APi  —  =  APlFui 


^<«2Ap2>  1  ^<u22rAp2>  + 


- +  - - -  w_  _u .  +  ,  _u_  - 

dt  r  dr  3/2  1  3/2  2 


(  16c) 


-Ap2fv2  *  AP2^_  "  Ap2Fu2 


W5 


ywmjt«JTKrMJ»KJnUltJntJ^Tni 


n.m  uivnuna  »»»  trm  .arm: 


d<  »0aPo)  i  a<  «,>v»is>  +  _  “  „  .  ^Q-gZg 

- _  +  _  .  +  w1/2vQ  W1/2V1 


<  17  »> 


+AP0fu0  *  Ap0Fv0 


d<v.Ap,,>  1  dtu.v.Ap.r)  +  -  + 

- i-S-  4 - -1-----  +  "3,2*1  -  w3/2v2  -  w1A 


1/2V0 


<  17b> 


AP1u,v1 

+  W1/2V1  +  ~ ~~  +  Aplful  *  APlFvl 


d< v2AP2>  i  a<  u2v2Ap2r> 
at  r  dr 


W3/2V1  +  *3/2v2  * 


AP2u2v2 


(  17c  > 


+aP2£u2  88  Ap2Fv2 


whtrt:  w  *  ff© 

.  1  -  1 
*#  =  — <  |w  I  +  W)  r  w  =  ~<  Iwl  -  w> 

2  2 


(  18  > 


(  19  > 


APj  “  p j  — 1/2  ~  Pj  +  l/2  *nd  P*  *  P-l/2  C"0> 

Not*  that  0  is  aaaum*d  to  b*  **ro  at  th*  aurfaca  and  at  th* 
uppar  boundary. 

Intagrating  <4>  acroaa  aach  layar  givaa: 


a<Ap0)  1  a<rUgApg)  ^ 

at  r  dr 


wl/2  *  0 


(21) 


a  •  «  "  *  * 


( .  f«*.  «*_ 


d<Ap 


1>  1  dtrUjAp,^ 

r  dr 


+  w3/2  “  wl/2  =  0 


(22) 


d<Ap2>  1  d<ru2Ap2> 

dt  r  dr 


“  w3/2  =  0 


(  23> 


where: 


w  =  w  -  w 


Ones  the  squat  ions  are  vertically  discretized*  it  ie 


more  convenient  to  work  with  the  momentum  equations  in 


advective  form. 


Multiplying  (21)  by  un  and  subtracting 


from  < 16a>*  and  multiplying  <21>  by  vn  and  subtracting 


from  < 17a)  gives: 


So  +  u  fHg  +  wl/2<u0  -  V  !0_  _  £v  +  ^0 

dt  0  dr  Apn  r  °  dr 


Fu0  <24> 


0  +  ,  3V0  .  Wl/2<v0  ~  vl)  ,  U0V0  ,  . 

4  uQ—  4 - - - 4 - 4  fu0 

dt  dr  Apn  r 


Fv0  <25> 


The  results  are  similar  for  layers  1  and  2. 


Provided  that  w+,  w  *  and  F  are  known*  then  <!!>-<  13), 


<21>-<25>  form  a  closed  system  if  the  p^  quantities  can  be 


related  to  M^. 


Def ining 


A©,  -  e0 


A0l/2  =  ei  e0 


(26) 


A03/2  =  02  “  °1 


v.v ’a.v.v.V'.v.'V.v.v.vv.a  -'././.Av.v.y.v.v,/ 


than#  from  <11>— <13># 


M0  ■  A®.n. 


B1  *  A®.  .  +  Aei/2nl/2 


«2  °  A®«n«  +  A®l/2ni/2  +  A®3/2n3/2 


From  (20) 


Ap0  *  P.  -  P1/2 


AP1  *  Pl/2  p3/2 


Ap2  =  P3/2  ~  PT 


and  so# 


p3/2  *  Ap2  +  PT 


pl/2  *  Ap2  +  Apl  +  PT 
pb  =  AP2  +  Apl  +Ap0  +  P1 


Now#  from  tha  chain  rula# 


an,  an,  dp, 

dr  dp  dr 

dn. an..,  dp 


1/2  atfl/2 


dp  dr 


dn3/2  _  dI13/2  ^£3/2 


dp  dr 


FljWWfl 


( 34  cont. > 


d<Ap_> 


<YS  +  Tl/2  +  73/2>  “ 


hr 


It  is  convenient  to  divide  7 ^  into  a  mean  part  y  ^ 
calculated  from  constant  pressure,  P^,  and  a  perturbation 


part  7j*  where 


m  yi  +  V 


<  35) 


Now*  define 


AP0  = 


—  €0h0 


Ap,  = 


:  *xhi 

S 


(36) 


Ap, 


=  Z  *2h2 

y _ 


m 

o 

II 

7 

s 

/  7 

8 

=  1 

^0  = 

7 

s 

/  7,* 

•l  = 

\ 

+  ^l/Z* 

7 

s 

/<  7  * 

B 

+ 

7  • 

1/2 

>  ( 37) 

*2  = 

”*8 

«yB 

+  ^1/2  +  73/2 

II 

<aa 

>s 

/<  7g’ 

+ 

7  1 

1/2 

*  Y3/2> 

where 

‘i 

represent  static 

stabi 1 ity 

factors 

of 

the  mean 

atmosphere*  the  inverse  of  /3j  represent  static  stability 
change  due  to  deviations  from  the  basic  state  of  the 
atmosphere*  and  the  factors  are  proportional  to  h^  the 


thickness  of  the  layers. 


Using  < 35)  and  <36),  (34)  can  be  written  as 


d<  11q  +  e^h^  +  «2h2*  ®  d<  +  *2^1  +  t2ix2) 


(38) 


d<  h^  +  h^  +  <«2/«^>h2>  g  dh^  g  £<  e^h^  +  *2h2> 


0O  dr  /3X  dr 


d<hQ  +  hj^  +  h2>  g  dhQ  g  d<  e^h^)  g  d(  f2h23 


0O  dr  dr  32  dr 


Since  t ^  and  /3^  are  nondimensional#  h^  have  the  dimensions 


of  length. 


Define  the  heating  term*  Q*  as 


7swj-l/2 


'j-1/2 


j-1*  2 


<  39> 


g  d<  hQ  +  ‘jhj^  +  «2h2> 


g  dhQ  g  d<  fjh^  +  *2h2> 


<  40> 


0n  dr  fi. 


g  dhQ  g  d<  «1h1>  g  d<  <2h23 


0q  dr  Px  dr 


*2 


and  the  geopotential*  4> *  as 


4>n  =  g  <  hQ  +  «1h1  +  <2h2> 


<t> ,  =  g  <hQ  +  h  ^  +  <«,/«,)  h,  > 


2  1  2‘ 


(41) 


wflBvwvvwy 


(41  cont .  > 


*2  =  g  <h0  +  hx  +  h2> 


than  <24>  can  be  written  as 


3U0  3“o  °l/2  to  3*o 

—  +  «o— +  <  uo  -  v - fvo +  — 

St  dr  hQ  r  Sr 


'  -%  *  Fu0 


a“l  SUj  83/5  °3/2,  V1 

+  «!-—  +  ~~ <“l  -  u2>  *  7~ <  "1  -  V - fv1 

St  dr  *1^1  *1^1  r 


+  sT m  +  Ful  <42> 


Su2  ,  °3/2. 


— -  +  u_— -  +  ---=<  u_  -  u.)  -  ---  -  fv_  +  — - 
dr  Sr  «0h,  r  Sr 


-*2  +  Fu2 


and  (25>  can  be  written  as 


3v0  .  avo  A  Qj/2, 

St  °Sr  hQ 


^  «. 
dt  Sr 


U0V0 

Vi>  +  --  +  fuQ  =  FvQ 


,  *  °3'2,  .  .  ulvl 

V  +  — --<v1  -  v0>  + - 

*lhl  r 


^v2  ®3/2 

+  u2~r  +  : ( v2 
dr  *2h2 


+  £u,  = 


U2V2 

Vi>  +  -±-±  +  £x,2 
r 


(  43) 


•c*v 


iy»V 


fv2  + 


u2  d<ru2>  u2  +  v2 

2  r  dr  r 


(  45  cont . > 


.  °3/2<ul  “  V 
+  -  +  F 


«2<H2  +  V* 


«0 

a<rV 

^  q1/2<v1  *  V 

r 

dr 

< H0  +  hQ’ > 

+ 

^1 

d(  rv^> 

„  °1/2<V0  "  V 

r 

dr 

*1<H1  +  h1*  > 

°3/2<v2  "  vl> 


<  46) 


+  h1*  > 


^2  d(rv2)  +  °3/2<vl  ~  V  +  p 


r  dr 


e2<H2  +  h2’> 


d<  ruQ> 

dh0’ 

d<  ruQ> 

_  k  1 _ _ 

+ 

.  n 

r  dr 

*  _uo“~~ 

dr 

u 

<T> 

U 

o 

m 

t 

i 

°l/2 

d<  ru. ) 

dh  • 

<3(  ru. > 

+ 

Qi 

d<  ru.) 


d<rV  .  °3/2  °3/2 


W 


M' 


where  and  <f>^  art  dafinad  in  aquations  (40>  and  <41># 
raspact ively.  Nota  that  tha  aquations  for  a  rotating 
ax isymmatr ic  fluid  in  isentropic#  cylindrical  coordinates 
ara  similar  to  tha  aquations  presented  by  Ooyama  (1969) 
for  an  atmosphere  with  three  incompressible  fluid  layers. 
Tha  only  differences  are  that  the  terms  do  not  appear 
in  the  incompressible  fluid  layer  model#  and  Ooyama* s 
model  uses  the  balance  approximation. 

In  this  3-layer  model  structure#  0Q,  0^,  02#  and 

are  specified  constants.  It  is  also  necessary  to  specify 
the  basic  state  pressures  pB#  P1/2*  and  p3/2*  0ncm  these 
values  are  specified#  the  constants  HQ#  H1#  H2#  *x#  and  <2 

can  be  calculated  as  follows: 


m  •• 

lk)K  •••- 


<  48) 


P3/2\  p0  / 


<  02  -  0X) 


=  1 


*1  s  \  n\  +  w 

* *>  =  t ^  t  7,  .«  ■+■  7«  i«) 


(  49) 


‘  P»  -  fl/2*  \ 


<Pl/2  ~  *3/2*  7s 


(  50) 


H  <P3/2-PT>7. 

H2  " 

9 

Tha  0^  tarns  in  (40)  ara  calculatad  from  h* q#  h' ^  and 
h' 2  and  tha  following  aquations: 


*3/2  *  *T  +  ‘•I£<H2  +  V>  *  *3/2 


®‘2 

+  — =<h  ’  > 

y_  2 


*1/2  "  *3/2  +  V*  “  *1/2  +  flhl’  +  ‘2V 


<  51) 


*s  s  *i/2  +  r^Ho  +  V>  *  p.  +  r< V  +  *ihi’  +  ‘2h2f> 


Rd/p. 


T  '  =  T_  -  T_ 


P.\*o 


ld  1*1/2 


*1/2  \  P0 


*1/2  \*0 


Tl/2'  “  Tl/2  ”  7l/2  <52> 


T  •  ■  T  —  T 

3/2  3/2  t3/2 


In  order  to  close  the  system  of  equations#  <40>-<41> 


end  <45>-<53>#  the  F  terms  snd  the  dlebatic  terms#  ^1/2 
and  Q3/2'  ■ua*  he  parameterized  in  terms  of  other 
dependent  variables. 

To  parameterize  F#  the  surface  drag#  horizontal  eddy 
diffusion  of  momentum#  and  the  vertical  shear  between  the 
layer  Interfaces  are  considered.  The  F  terms  in  equations 
<  54)  ere  written  as  ( DeMar ia  and  Schubert#  1984): 


F0  * 


CD -  2  1  - 

f  iv0ivQ  4.  x2<vz  -  -2>v0 

Ho  r 


-  M 


(W 

<H0  +  h0*> 


(54) 


Vv  -  Vvi 


M<VV 


•f  - 


VV 


*1<H1+h1,>  *1<H1+h1’  > 


F2  = 


2  1  _  "<W 

X  (VZ  -  -  >V_  ♦ - 

WV* 


<  A) 


<  B) 


<C> 


where  IVq  |  represents  the  magnitude  of  the  boundary  layer 
horizontal  velocity.  Term  < A>  represents  the  surface  drag 


calculated  from  tha  bulk  aarodynamic  formula.  Am  in 
DaHaria  and  Bchubart  < 1984 >  and  Ooyama  <1969)*  tha  drag 
coafficiant  CD  is  assumed  to  have  constant  value  of  0.0015. 
Terms  < B)  represent  horizontal  eddy  diffusion  of  momentum. 
The  horizontal  eddy  diffusion  coefficient  was  assumed  to 
have  a  constant  value  of  103  m2/m ,  M  ln  D.Maria  and 
Schubert  <1984).  Terms  < C)  result  from  vertical  diffusion 
caused  by  velocity  shear  across  the  interface  between 
layers  0  and  1  and  layers  1  and  2. 

To  parameterize  Q^/2'  assume  that  the  pressure 
thickness  of  the  boundary  layer  < layer  0)  remains  constant 
with  time  so  that  dh/dt  *  0.  Then  from  <47>* 

-  +  d<ru_> 

°l/2  “  °l/2  *  <**0  +  h0,?  Z  <57> 

r  or 

+ 

Since  *nd  Qj/2  mrm  allfaF*  positive*  then 

*0  <*  0 
«0  >  0 

<  58) 

«0  <  0 

60  >=  0 

1  d<ruQ) 

where  *  -  -  *  horizontal  divergence  in  layer  0. 

r  dr 

+  - 

Parameterizing  Q^2  and  ®i/2 


using  <58)  indicates 


that  all  of  tha  uai  which  converges  in  layer  0  is 
transported  to  layer  1  by  diabatic  heating  and  that  mass 
divergence  in  layer  0  is  replaced  by  mss  from  layer  1  due 
to  diabatic  cooling. 

To  parameter  ice  Q3/2#  assume  that  this  term  represents 
diabatic  heating  due  to  cumulus  convection.  Since 
sustained  cumulus  convection  requires  boundary  layer 
moisture  convergence#  then  assuM  that 
+  + 

°3/2  '  «l/2  <59> 

Equation  (59)  indicates  that  for  each  unit  mss  of  air 
transported  from  layer  0  to  layer  1#  y  units  of  mass  are 
transported  to  layer  2.  OoyaM  (  1969)  noted  that  y  is  a 
measure  of  instability  of  deep  convection  that  includes 
the  stabilizing  effect  of  the  warm  core  and  the  variation 
of  equivalent  potential  temperature  in  the  planetary 
boundary  layer  < layer  0).  The  moist  static  energy  <S*2> 
of  the  saturated  cloud  air  which  reaches  layer  2  is  given 
by: 


where 


_  *  _  8o  +  ~  1>si 

B2  „ 


S  = 


S  = 


s  = 
* 

s  = 


CpT  +  gz  +  Lq 


M  +  Lq 


ft  ft 

CpT  +•  gz  +  Lq  =  M  +  Lq 

moist  static  energy  of  unsaturated  air 

moist  static  energy  of  saturated  air 


(60) 

(61) 
(  62) 


L  =  latent  heat  of  vaporization 


q  =  mixing  ratio  of  unsaturatad  air 
* 

q  =  mixing  ratio  of  aaturatad  air  < saturated 
mixing  ratio) 

=  vertical  stability  parameter 
Solving  (60)  for  »?  gives: 


The  Montgomery  stream  functions  in  <61>-<62>  can  be 
calculated  using  the  following  derived  from  <7>  and  (27): 


M 


1 


M 


2 


(64) 


The  Montgomery  stream  functions  which  appear  in  the 
definitions  of  and  Sj  do  not  represent  the  values  of  M 
in  the  model  layers*  but  rather  represent  M  at  the  pressure 
levels  of  entrainment  PRE  and  detrainment  PRD.  These  are 
linearly  interpolated  from  the  M  values  of  the  model  layers 
assuming  that  these  are  representative  of  the  value  of  M  at 
the  pressure  of  the  mid-points  of  each  layer  Pf  and  Pr, 
which  are  given  by: 


25 


(  65) 


Fl  *  <Pl/2  +  p3/2>  1  2 


p0  =  ' p3/2  +  PT 


+  P*>  /  2 


Let  the  Montgomery  stream  functions  at  p^  and  p^  be  M^  and 
M2*  end  the  Montgomery  stream  functions  at  PEE  and  PRO  be 
Mg  and  Mp.  Then  Mg  and  Mg  are  linearly  interpolated  from  pL 


and  pn  by  the  following: 


<M  -  M.  > 

«E  s  M1  +  ~ - ^  PL  "  PRE> 


<  PL  "  Pu> 


(  67) 


<M?  -  M. > 

md  =  Ma  +  — - - --<  pL  -  PRO) 


< PL  "  P0> 


Then  the  moist  static  energies  are  calculated  UBing: 


S0  s  H0  +  L*r 


S1  -  Mg  +  Lq. 


<  68) 


S2*  *  M0  +  L^2* 


The  moisture  values  and  q  2  could  be  calculated  at  each 

A 

time  step.  However*  the  variations  of  Lq^  and  Lq  2  are  not 
large  compared  to  changes  in  and  S*2 *  respectively*  and 

H 

so  q ^  and  q  2  will  be  held  constant.  Variations  of  the 


boundary  layer  moisture*  qn*  however*  can  be  large.  For 


Ls  reason*  qQ  is  predicted  using  the  following 


conservation  equation: 


+  “or 

dt  dr 


<Hq  +  h0’> 


C„|V^| 


<  Hq  +  hQ*  > 


where  q#  is  ths  saturation  mixing  ratio  at  T>b  and  p^;  T 
bsing  ths  sea  surface  temperature*  and  pB  the  surface 
pressure.  Although  Mg  is  a  function  of  p^*  it  can  be  seen 
from  <64>  that  Mg  will  decrease  as  pB  decreases. 
Observations  show  that  as  a  storm  intensifies*  the 
temperature  at  the  surface  remains  approximately  constant 
despite  the  decrease  in  the  surface  pressure.  This  is  due 
to  the  sensible  heat  flux  from  the  ocean  surface.  Since  it 
is  assumed  that  9  =  0  at  s  =  0*  the  surface  sensible  heat 
flux  cannot  be  included  in  the  isentropic  coordinate 
model.  The  use  of  the  predicted  Mg  might  then  lead  to 
artificially  low  values  of  Sg.  To  compensate  for  the  lack 
of  surface  sensible  heat  flux  in  the  model*  the  value  of 
Mg  of  the  reference  state  ,  Sg,  will  be  used  in  the 
calculation  of  Sn  rather  than  the  predicted  value.  Then 


=  Mn  +  Lgf 


where 


Vo1 


(  70) 


(71) 


The  remaining  term  to  be  parameterized  is  Q3/2* 
This  term  represents  radiative  cooling.  If  FgLR  and  FgLp 
represent  typical  radiative  cooling  rates  in  clear  and 

cloudy  regions*  respectively*  then 


zero-  or  first-order  modified  Bessel  function  of  the 
first  kind  for  the  u*  v*  and  h'  variables  or  zero-  or 
Bessel  functions  for  the  qQ  variable.  The  type  of  basis 
function  used  is  dependent  upon  the  lateral  boundary 
conditions  <DeMaria*  1983).  In  this  case*  for  determining 
the  basis  functions  for  u*  v,  and  h' *  geostrophic  balance 
is  assumed  at  r  =  a*  the  outer  boundary  of  the  model*  and 
u=v=0atr=0*  the  center  of  the  storm.  The  lateral 
boundary  conditions  for  are  that  the  radial  gradients 
of  gQ  are  zero  at  r  =  0  and  a. 

Once  the  horizontal  transform  has  been  applied*  the 
governing  equations  are  solved  using  second-order 
Adams-Bashforth  time  differencing  with  an  initial  forward 
time  step.  In  the  simulations  presented  in  the  following 
sections*  the  horizontal  series  are  truncated  at  N  =  32 
with  a  =  800  km.  The  nonlinear  terms  of  the  governing 
equations  were  calculated  using  the  transform  method 
< Orszag*  1970;  El iasen  et  al*  1970)  using  3N  +5  evenly 
spaced  grid  points.  Therefore*  the  spacing  of  the 
transform  grid  is  8  km. 

DeMaria  and  Schubert  ( 1984)  and  Machenhauer  ( 1979) 
described  the  phenomenon  of  spectral  blocking  and  the  use 
of  linear  diffusion  terms  to  mitigate  the  errors  associated 
with  it.  Spectral  blocking  is  a  result  of  the  higher 
wavenumbers  of  the  truncated  spectral  series  increasing  in 
energy  caused  by  the  steepening  gradients  of  the  dependent 


1 


variables  near  the  storm  center.  In  this  study  it  was 

4 

found  that  linear  diffusion  terms  of  the  form  v 
effectively  inhibited  spectral  blocking.  Terms  of  this 

form  with  =  1.671  x  10 ^  m^ /s  were  added  to  all  of  the 

prognostic  governing  equations.  In  a  simple  one¬ 

dimensional  diffusion  problem#  the  amplitude  of  highest 
mode  in  a  Fourier  series  with  32  modes  on  an  800  km 

periodic  domain  decreases  by  a  factor  of  e  in  4  minutes 
using  this  form  of  diffusion. 

2.5  Initial  Conditions 

The  initial  mass  and  wind  fields  for  all  simulations 
run  in  the  study  are  given  by: 
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where  the  maximum  tangential  wind#  v_,  equal  to  15.0  m/s  at 

ill 

a  radius#  r  #  of  100  km.  This  initial  condition  is  an 

III 

axlsymmetric  vortex  in  gradient  balance  with  no  initial 
radial  circulation.  The  values  of  y 

30 


m 


and  r_  were -chosen 
in 


after  preliminary  model  runs  showed  a  realistic  and  cost 
effective  rate  of  intensification  (reaching  maximum 
intensity  in  approximately  4  or  5  days)  and  a  smooth 
development  in  the  radius  of  maximum  tangential  winds. 

To  run  the  model,  it  is  necessary  to  specify  the 
following  parameters: 

PT'  p3/2'  pl/2'  Ps'  ®2'  ®1'  90'  and  Tss 5 
and  initial  conditions  must  be  specified  for 

u0'  ul'  u2'  v0'  vl'  v2'  h0'  hl'  h2'  *nd  q0- 

Table  1  shows  the  numerical  constants  used  in  the 
numerical  simulations.  Average  values,  such  as  for  g  and 
6,  were  derived  from  average  soundings  of  typical  tropical 


Table  1:  Values  of  parameters  In  the  standard  numerical 


model  simulation. 

NUMBER  OF  NORMAL  MODES 

N 

32 

RADIUS  OF  STUDY 

AX 

800  km 

TIME  STEP 

DT 

60  s 

LENGTH  OF  SIMULATION 

5  days 

SECOND-ORDER  DIFFUSION 

DA2 

1.0*  103 

FOURTH-ORDER  DIFFUSION 

DA4 

1.  67 1*10+13 

P0 

1000  mb 

P. 

1010  mb 

Pl/2 

910  mb 

P3/2 

500  mb 

P5/2=PT 

150  mb 

ENTRAINMENT  PRESSURE  LEVEL 

PRE 

700  mb 

DETRAIN MEN T  PRESSURE  LEVEL 

PRD 

300  mb 

eo 

300°X 

©1 

315°X 

^2 

345°X 

LATITUDE 

20°N 

°P 

1000  J/kg°X 

X 

0. 2859 

9 

9.8 

M 

5xl0“4 

CE 

0.0015 

CD 

0.0015 

*0 


20  gm/kg 


Table  1  (continued) 


9^  6  gm/kg 

q£  0.8  gm/kg 

Tsb  28°C 

Lv  2.5xl06  J/kg 


The  standard  modal  simulation#  for  which  th®  ®f facts  of 

varying  sea  surface  temperature#  mid-level  moisture#  and 

latitude  were  compared#  had  a  sea  surface  temperature  of 

28°C#  mid-level  moisture  of  6  g/kg#  and  a  latitude  of 
o 

20  N.  The  model  simulation  represents  5  days#  as  did  all 
subsequent  runs.  The  amount  of  CPU  time  used  for  each  run 
was  roughly  3.5  hours  on  a  VAX  11/750. 

Figures  2  and  3  represent  the  radial  and  tangential 
wind  profiles  for  layers  0#  1#  and  2  at  t  =  0#  48#  and  96 

hours.  Winds  in  layer  0  rapidly  increase  the  first  two 
days#  level  off  during  days  3  and  4#  and  decay  during  day 
5.  This  corresponds  to  the  daily  minimum  surface  pressure 
trends  (fig.  4>.  Tangential  wind  profiles  in  layer  1  are 
very  similar  to  layer  0  except  slightly  slower.  Winds  in 
layer  2  are  cyclonic  near  the  center  but  become 
anticyclonic  at  distances  roughly  200  to  400  km  from  the 
storm  center.  Radial  winds  in  layer  0  are  primarily 
convergent  and  divergent  in  layer  2  (fig.  5).  The  zones  of 
maximum  convergence  in  layer  0  and  maximum  divergence  in 
layer  2  are  between  50  to  100  km  from  the  storm  center. 

Figure  6  shows  the  profiles  of  n#  the  vertical 
stability  factor.  Decreasing  values  near  the  storm 

center  and  at  r  -  100  km  are  due  to  subsidence#  which  dries 
the  boundary  layer  and  reduces  SQ.  High  values  near  r  =  75 


RADIUS  (km) 


Figure  2  Profiles  of  tengentisl  winds  in  leyers  0,  1,  and 
2  et  t  ■  0  (solid),  46  (dash)#  and  96  (dash-dot) 
hours. 


m/s 


RADIUS  (km) 


Figura  3  Prof  lias  of  radial  winds  in  layars  0*  1#  and  2 
at  t  *  0  (solid)*  48  (dash)*  and  96  (dash-dot) 
hours. 


Figure  5  Profiles  of  divergence  in  layers  0  and  2  at  t 
0  (solid)#  48  (dash)#  and  96  (dash-dot)  hours. 
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km  act  a  raault  of  increased  avaporation  as  tha  dry# 
subsiding  air  at  r  =  100  km  movas  laterally  toward  the 
storm  canter.  The  overall  gradual  decrease  of  n  in  the 
central  region  of  the  storm  is  caused  by  the  formation  of 


the  warm  core#  which  inert 


®2  ~S^  with  time.  The  outer 


regions  of  the  storm  show  decreasing  ij  values  with  time 
caused  by  large  scale  subsidence  resulting  from  low-level 
divergence#  which  tends  to  reduce  the  moisture  in  the 
boundary  layer. 

The  radii  of  maximum  radial  and  tangential  winds 
increase  for  the  first  12  hours#  then  gradual  decrease 
until  day  2.5#  followed  by  slow  expansion  until  the  end  of 
the  simulation  (figs.  7  &  6).  The  radii  of  30-  and 
50-knot  winds  were  used  by  Weatherford  ( 1985)  to  describe 
the  aerial  size  of  storms.  In  metric  units#  these  represent 
the  radii  of  15  and  25  m/s  winds#  respectively.  The  radii 
of  30-  and  50-knot  winds  in  layers  0  and  1  rapidly 
increase  during  the  first  hours  after  the  storm  attains 
winds  of  30-  and  50-knot  winds  (fig.  10).  Following  the 
rapid  expansion#  the  size  of  the  storm  increases  by  steady# 
near-linear  growth.  The  radius  of  30-knot  winds  increases 
faster  than  the  radius  of  50-knot  winds. 

The  basic  state  temperature  profile  used  in  the 
previous  simulation  represented  a  conditionally  unstable 
atmosphere  and  the  boundary  layer  was  fairly  moist 
initially  (relative  humidity  «  82%).  Since  a  hurricane  is 


v max  (km) 


a  fairly  rare  event#  a  tropical  cyclone  model  ehould  be 
able  to  simulate  nondevelopment  as  well  as  development 


(Emanuel#  1986).  To  test  this  case#  the  basic  state 
potential  temperature  structure  was  changed  to  represent  an 
atmosphere  which  was  exactly  moist  adiabatic  < -  300°K# 

6 2  -  316°K#  and  #2  =  350°K) #  and  the  initial  boundary 

layer  moisture  was  reduced  to  qQ  =  16  gm/kg  (relative 
humidity  =  67%). 

Simulations  with  maximum  tangential  wind  speeds  of  5 
m/s  produced  no  intensification  in  5  days#  whereas  initial 
conditions  with  winds  of  15  m/s  produced  a  strong  tropical 
storm.  The  15  m/s  case  was  still  intensifying  after  5  days 
of  simulation#  and  reached  a  maximum  intensity  of  991  mb  in 
6  1/2  days.  Using  the  conditionally  unstable  temperature 
and  moisture  profile#  the  5  m/s  case  intensified  to  1008 
mb.  These  results  are  qualitatively  similar  to  Rotunno  and 
Emanuel's  (1987)  simulations  which  used  an  axisymmetric 
nonhydrostatic  numerical  model  which  explicitly  resolved 
the  cumulus  convection. 

3.2  Comparison  of  Cumulus  Parameterization  Schemes 

The  major  difference  between  the  study  model  and 

Ooyama's  model  (1969)  is  the  cumulus  parameterization 

scheme.  Both  schemes  involve  the  calculation  of  the 

vertical  stability  factor#  *7  Ooyama  assumed  c^t/6  is 

* 

approximately  constant  so  SQ#  S^#  and  S2  in  the  definition 
of  7  can  be  replaced  by  0#Q  ,  0#1#  and  0#2*  respectively. 
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is  predicted  using  a  conservation  equation 
that  used  to  calculate  dqQ/dt#  0^  is  assumed 

A 

constant#  and  6^  is  diagnosed  from  the  thickness  of  layer 

I  2  using  an  analogy  with  a  compressible  fluid  and  an 

\ 

|  empirical  relation  between  upper  level  0^*  changes  and 

» 

!  mid-level  changes.  The  parameterization  scheme  in  this 

study  allowed  for  more  specific  control  of  measurable 
|  parameters  in  the  calculation  of  >?#  such  as  sea  surface 

temperature  and  moisture  content  of  the  layers. 

v  Keeping  the  initial  conditions  similar  (table  2>#  the 

* 

l 

I  two  cumulus  parameter izat ions  were  compared.  The  results 

J  from  each  parameterization  scheme  were  qualitatively 

similar.  The  major  differences  were  that  simulations  using 
[  Ooyama's  parameterization  intensified  faster  than  (reaching 

maximum  intensity  almost  two  days  before)  but  not  to  the 
extent  of  (14  mb  less)  the  study  model  (fig.  10).  In 
|  addition#  Ooyama's  model  produced  a  slightly  larger  storm# 

based  on  the  radius  of  30-knot  surface  winds  (R-30). 

3. 3  Sensitivity  Studies 
I 

The  effects  of  individual  parameters  in  tropical 

l 

!  cyclone  development  were  studied  by  systematically  varying 

[  one  parameter  during  model  runs  while  keeping  all  other 

! 

parameters  constant.  Parameters  that  were  varied  are  sea 

i 

surface  temperature#  mid-level  moisture#  and  latitude  of 
the  storm.  The  sea  surface  temperature  ranged  from  24°  to 

!  o 

30  C#  mid-level  moisture  from  2  to 


Then  0__ 

BO 

similar  to 


10  gm/kg  (relative 


Table  2:  Vertical  equivalent  potential  temperature 

profiles  used  in  testing  of  Ooyama’s  (1969)  model 
and  the  study  model. 
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359°K 

Figure  10  Comparison  of  results  from  the  study  numerical 
model  <  solid)  and  Ooyama's  (  1969)  numerical 
model  (dashed):  (a)  minimum  surface  pressure 

(SLP)  with  time*  ( b)  radius  of  30-  and  50-knot 
surface  winds  ( R-30  &  R-50)  with  time*  and  (c> 

maximum  radial  versus  tangential  winds  (both 
model  results  are  identical). 


humidity  from  16  to  82%) #  and  latitude  from  5°  to  40°N. 
The  effect  of  seasonality  was  studied  by  using  latitudinal 
monthly  average  sea  surface  temperatures  for  July* 
September#  and  November  in  the  western  Pacific  calculated 
from  Crutcher  and  Davis  <1969). 

The  basic  characteristics  used  to  describe  the  state 
of  the  tropical  cyclone  were  the  minimum  surface  pressure# 
rate  of  intensification  (hourly  difference  of  the  central 
surface  pressure)#  maximum  wind  speed*  radius  of  maximum 
wind  speed*  and  outer  radius  of  30-)cnot  and  50-knot 
surface  winds  (used  as  a  measure  of  cyclone  size  in 
Weatherford#  1985).  Individual  model  runs  were  compared 
for  given  time  intervals  from  initiation  and  at  two  stages 
in  the  cyclone  development:  the  period  of  minimum  surface 
pressure  (maximum  intensification)  and  the  period  of 
maximum  hourly  change  of  the  central  surface  pressure. 

3.3.1  Effect  of  Sea  Surface  Temperature 

The  minimum  surface  pressure  (maximum  intensity) 
observed  during  the  5-day  simulation  decreases  as  the  sea 
surface  temperature  increases  (fig.  11).  The  effect  of 
sea  surface  temperature  on  the  maximum  cyclone  intensity 
is  dominant  compared  to  the  effects  of  mid-level  moisture 
and  latitude.  The  maximum  rate  of  intensification 
increases  nonl inear ly  (fig.  12)*  while  the  time  it  takes 
the  storm  to  reach  maximum  intensification  from  initiation 


deer 


linearly  with  increasing  sea  surface  temperature 


SLP  (mb) 


Figure  11 
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Minimum  surface  pressure  (SLP)  at  maximum 
intensity  with  respect  to  sea  surface 
temperature  (solid)*  mid-level  moisture  (dash- 
dot)*  and  latitude  (dash). 
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(fig.  13>.  In  addition#  storms  that  have  developed  over 
warm  waters  appear  to  deintensify  more  quickly  than  storms 
over  cooler  waters  (fig.  14a). 

Figure  14b  shows  the  rate  of  growth  in  storm  size. 
Starting  from  initiation#  storms  developing  over  warmer 
waters  grow  larger  and  more  rapidly  than  storms  developing 
over  cooler  waters.  However#  because  the  storms  intensify 
more  rapidly  in  warmer  waters  the  storms  are  roughly  the 
same  size  at  the  period  of  minimum  sea  surface  pressure. 
Note  that  storms  continue  growing  in  size  during 
deintensification  (filling).  The  only  case  in  which  the 
size  of  the  storm  decreased  during  the  five  day  simulations 
was  when  sea  surface  temperature  was  decreased  during  the 
time  integration.  The  effects  of  changing  6ea  surface 
temperature  were  estimated  with  this  model  by  changing  the 
sea  surface  temperature  instantly  for  the  entire  storm 
domain.  The  effects  of  changing  sea  surface  temperature 
drastically  and  gradually  were  essentially  the  same:  the 
storms  deintens if ied  and  decreased  in  size  when  the  sea 
surface  temperature  was  reduced. 

3.3.2  Effect  of  Mid-Level  Moisture 

The  effect  of  mid-level  moisture  is  difficult  to  model 

since  q^  ig  held  constant  throughout  the  integration. 
However#  it  is  worthwhile  to  determine  if  the  model 


sensitivity  is  qualitatively  reasonable.  Although  the 
mid-level  moisture  affects  the  cyclone  development  through 


Figure  14  Effect  of  sea  surface  temperature  < BST>  on  the 
<a>  minimum  surface  pressure  ( SLF)  with  time  and 
on  the  < b>  radius  of  30-knot  winds  in  layer  0 
<  R-30)  with  time. 


the  variable#  v,  as  does  the  sea  surface  temperature#  the 
influence  of  mid-level  moisture  on  tropical  cyclone 
development  is  not  identical  to  that  caused  by  sea  surface 
temperature.  Figure  11  shows  a  linear  relationship 
between  maximum  intensity  and  mid-level  moisture  whereas 
the  sea  surface  temperature  variations  produced  nonlinear 
variations  in  intensity.  Although  changing  mid-level 
moisture  produces  smaller  changes  in  maximum  cyclone 
intensity#  changing  mid-level  moisture  creates  large 
changes  in  the  maximum  rate  of  intensification  compared  to 
sea  surface  temperature  (fig.  12). 

In  this  model#  mid-level  moisture  and  sea  surface 
temperature  generally  influence  tropical  cyclone 
development  similarly.  Storms  intensify  more  rapidly  under 
favorable  conditions#  and  reach  a  lower  minimum  surface 
pressure.  Despite  growing  larger  and  stronger  more  rapidly# 
at  the  time  of  maximum  intensity  the  storms  are  roughly 
the  same  size  no  matter  how  moist  the  mid-levels  are  (fig. 
15).  In  addition#  storms  with  moist  mid-levels 
deintensify  faster  than  storms  with  drier  mid-levels# 
similar  to  storms  with  a  warmer  sea  surface  temperature. 

3.3.3  Effects  of  Latitude 

Latitude  affects  the  development  of  tropical  cyclones 
through  the  CorioliB  parameter#  f#  which  appears  in 
equations  (45)  and  (46);  whereas  the  effects  of  sea 


surface  temperature  and  mid-level  moisture  are  represented 


through  v »  the  vertical  stability  factor.  In  equations 

(45)  -  <  47 )  i\  is  included  through  the  Q+  an(» 

1/2  ana  w  3/2 

terms . 

With  all  other  factors  held  constant*  tropical 
cyclones  intensify  most  between  20°  and  30°N  (fig.  ID* 
although  storms  at  lower  latitudes  experience  greater 
rates  of  intensification  than  storms  at  higher  latitudes 
(fig.  12>.  Starting  from  initialization*  storms  reach 
maximum  intensity  in  roughly  3.5  days  independent  of 
latitude  (fig.  13)*  although  after  2  days*  the  cyclone  at 
10°N  did  not  intensify  significantly  (fig.  16).  Storms  at 
higher  latitudes  grow  larger  and  faster  than  storms  closer 
to  the  equator.  In  addition*  storms  at  maximum  intensity 
are  significantly  larger  toward  the  higher  latitudes  (fig. 


In  summary*  as  sea  surface  temperature  and  mid-level 
moisture  increase*  storms  intensify  more  rapidly  and 
intensify  to  a  greater  extent  (fig.  18) *  yet  the  storm  size 
is  the  same  when  at  maximum  intensity.  As  latitude 


i*  storms  intensify  less  rapidly*  generally 


intensify  to  a  greater  extent*  and  are  larger  when  at 
maximum  intensity. 


3.  3.  4 


rms.  hi*:  n  14)14  ■  i!!-a*ri?iCT-ri 


It  has  been  shown  that  sea  surface  temperature  affects 
the  intensity  and  size  of  tropical  cyclones*  but  are  the 
effects  similar  at  different  latitudes?  Figures  19  -  21 


Figure  16  Effect  of  latitude  on  the  <a>  minimum  surface 
pressure  < BLP)  with  time  and  on  the  ( b>  radius 
of  30-knot  winds  in  layer  0  (R-30)  with  time. 
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Figure  17  Radius  of  30-knot  winds  at  maximum  intensity 
with  respect  to  sea  surface  temperature  (solid)# 
mid-level  moisture  (dash-dot)#  and  latitude 
(  dash) . 
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Figure  19  Comparisons  of  trends  between  sea  surface 
temperature  < SST>  and  minimum  surface  pressure 
< SLP)  at  maximum  intensity  for  different 
lat itudes. 
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Figure  20  Comparisons  of  trends  between  sea  surface 

temperature  < 5ST>  and  the  rates  of  maximum 
intensification  at  different  latitudes. 
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Comparisons  of  trends  between  see  surface 
temperature  < 6ST>  and  radius  of  30-knot  winds  at 
maximum  intensity  <R-30>  at  different  latitudes. 


show  the  effects  of  sea  surface  temperature  on  storms  at 


different  latitudes.  Sea  surface  temperature  influences 
the  maximum  intensity  attained*  the  maximum  rate  of 
intensification*  and  storm  size  similarly  at  different 
latitudes*  yet  the  magnitude  of  the  effects  vary  with 
latitude.  The  nonlinear  effect  of  latitude  on  tropical 
cyclone  development  is  particularly  noticeable  when 
considering  maximum  intensity  and  maximum  rate  of 
intensification  attained  by  the  numerical  simulations 
(figs.  19  &  20).  In  the  case  of  maximum  intensity*  storms 
appear  to  become  strongest  between  20°  and  30°N  for  a 
given  sea  surface  temperature.  Storms  generally  display 
greater  intensification  rates  eguatorward*  but  this 
reaches  a  maximum  around  10°N.  At  5°N,  storms  do  not 
intensify  as  much  or  as  quickly.  The  weak  storms 
developing  at  S°N  may  be  a  result  of  poor  resolution  on  the 
8  km  grid.  Due  to  the  effect  of  sea  surface  temperature 
influencing  the  rate  of  intensification*  storm  size  at 
maximum  intensity  is  generally  independent  of  sea  surface 
temperature  and  highly  dependent  upon  the  latitude  (fig. 

21). 


3.3.5  Effect  of  Seasonality 

The  monthly-averaged  sea  surface  temperatures  for  the 
western  Pacific  (fig.  22)  were  used  to  estimate  the  effect 
of  seasonality  on  tropical  cyclones.  In  previous  sections* 
with  the  sea  surface  temperature  held  constant*  the  most 
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Figure  22  Latitudinal  profilaa  of  monthly-avaragad  aaa 
■ucfaca  tamparaturaa  < 8ST)  (from  Crutchar  and 
Davis#  1969). 


intense  storms  occurred  between  20°  an(j  30°N.  When  the 
latitudinal  dependence  of  sea  surface  temperature  is 
included#  the  most  intense  storms  occur  between  15°  and 
20°N  in  July  and  September#  and#  although  muoh  weaker# 
between  10°  and  20°N  in  November  (fig.  23>.  In  addition# 
storms  located  between  10°  and  15°N  intensify  the  fastest 
(fig.  24).  The  size  of  the  tropical  cyclone  at  maximum 
intensity  does  not  vary  significantly  during  the  season 
(fig.  25)  as  expected  from  the  lack  of  influence  of  sea 
surface  temperature  on  storm  size  at  this  stage  of 
development. 


The  numerical  model  in  this  study  predicts  trends  in 
the  size#  intensity#  and  rate  of  intensification  of 
tropical  cyclones  related  to  sea  surface  temperature  and 
latitude.  Climatology  data  of  tropical  cyclones  in  the 
western  Pacific  were  examined  to  see  if  these  predicted 
trends  in  size  and  intensity  caused  by  sea  surface 
temperature  and  latitude  are  realistic.  Data  gathered  by 
Air  Force  weather  reconnaissance  flights  in  the  western 
Pacific  from  1980  to  1982  and  reduced  by  Weatherford 
(1985)  was  used  to  study  the  effects  in  latitude  on 
tropical  cyclone  development.  Variables  that  were  used 
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Figure  23  Seasonal  changes  in  the  minimum  surface  pressure 
< SLP)  at  maximum  intensity  with  latitude. 
Dashed  line  represents  model  runs  with  sea 
surface  temperatures  of  28°C. 
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R-30  :  the  radius  (Jem)  at  which  the  35- 
knot  700  mb  wind  lay  on  average 
(approximately  the  30-knot  wind  at 
the  surface! 

LAT  :  latitude  of  cyclone  center 

INTL  :  past  intensity  change  in  meters. 

Difference  in  700  mb  minimum  heights  from 
past  12  hours  to  time  of  report 

SLP  :  minimum  observed  sea  level  pressure. 

The  data  were  analyzed  by  either  using  all 
observations  of  storms  regardless  the  stage  of  development 
or  using  data  of  storms  at  specific  stages  of  development. 
Data  that  included  multiple  observations  of  individual 
storms  were  categorized  by  intervals  of  latitude*  minimum 
observed  sea  level  pressure*  radius  of  30-kt  winds  and 
rate  of  change  of  latitude.  Using  the  categorized  data* 
mean  trends  were  compared  to  the  numerical  simulation 
predict  ions. 

When  the  data  of  individual  storms  was  complete  for  at 
least  several  days  of  observation*  specific  stages  of 
development  could  be  recognized.  Two  data  sets  were 
created  based  upon  the  following  stages: 

(1)  the  period  of  maximum  change  in  SLP  from  the 
observations  12  hours  before  and  after  that 
observation*  and 

(2)  the  period  of  minimum  SLP. 


When  the  observations  of  different  tropical  cyclones 

taken  throughout  their  development  were  compiled,  two 

general  trends  were  observed:  < 1)  storms  are  larger  at 

more  northerly  latitudes  (fig.  26),  and  (2)  storms  reach 

maximum  intensity  between  15°  an(j  30°N  (fig.  27).  Figure 

28  shows  the  effect  of  latitude  on  storms  of  different 

ranges  in  size  and  intensity.  Small-,  medium-,  and 

o 

large-size  storms  reached  maximum  intensities  between  15 
and  25°N.  The  model  results  for  the  case  when  sea  surface 
temperature  and  latitude  were  varied  agree  with  the 
climatological  results.  Without  considering  the  effect  of 
cooler  sea  surface  temperature  northward,  the  numerical 
model  simulations  reached  maximum  intensity  between  20°  and 
30°N,  north  of  the  climatology  results.  In  addition,  more 
intense  storms  are  larger  regardless  of  the  latitude.  When 
the  latitude  and  minimum  surface  pressure  were  considered 
(fig.  28b),  storm  size  increased  with  latitude  for  a  given 
intensity  category,  and  storm  size  increased  as  intensity 
increased  for  a  given  latitude  category. 

Because  numerical  simulations  and  climatology  Indicate 
a  strong  dependence  of  latitude  with  the  development  of 
tropical  cyclones,  the  influence  of  storm  movement  from  one 
latitude  to  another  was  studied  cl imatological ly.  The 
latitudinal  movement  of  storms  could  not  be  realistically 
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Figure  26 
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Climatology  of  tha  radius  of  30-knot  winds  ( R- 
30)  with  respect  to  latitude.  The  solid  line 
represents  the  mean*  and  the  dashed  line 
represents  one  standard  deviation  from  the  mean. 
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Figure  27  Climatology  of  the  minimum  surface  pressure 
(SLP)  with  respect  to  latitude.  The  solid  line 
represents  the  mean*  and  the  dashed  line 
represents  one  standard  deviation  from  the  mean. 


Figure  28  Climatology  of  < a)  the  minimum  surface  pressure 
< SLP)  with  respect  to  latitude  and  storm  size 
( R-30)  and  ( b>  the  radius  of  30-)cnot  winds  ( R- 
30)  with  respect  to  latitude  and  SLP. 


studied  using  the  axisymmetric  model  in  this  study.  The 
rate  of  latitudinal  movement  was  calculated  in  the 
following  ways:  (1)  the  difference  between  the  current 
latitude/  and  the  latitude  at  the  first  observation  divided 
by  time  between  the  two  observations  (called  the  average 
rate  of  latitudinal  movement)  and  (2)  the  difference 
between  the  current  latitude  and  the  latitude  12  hours 
previously  (called  current  rate  of  latitudinal  movement). 
The  average  and  current  rates  of  latitudinal  movement  were 
categorized.  The  effects  of  latitude  and  latitudinal 
movement  on  tropical  cyclone  size  and  intensity  were 
studied  (figs.  29  &  30).  Except  for  storms  moving 
southward/  storms  generally  reached  maximum  intensity 
between  15°  and  25°N  regardless  of  the  rate  of  poleward 
movement/  as  seen  in  the  numerical  simulations.  This 
trend  held  for  average  rate  of  latitudinal  movement  and  the 
current  rate  of  latitudinal  movement.  In  addition/  the 
size  of  the  tropical  cyclones  generally  increased  with 
latitude  regardless  of  the  rate  of  northward  movement 

4. 3  Climatology  of  Tropical  Cyclones  with  Respect  to 
Developmental  Stages 

In  order  to  study  differences  of  tropical  cyclones  at 
various  stages  of  development/  the  climatology  data  was 
separated  into  groups  based  on  maximum  rates  of 


intensification  or  maximum  intensity.  The  results  are  very 
similar  to  the  results  of  the  previous  climatology  section 


Figure  29  Climatology  of  (a)  minimum  surface  pressure 
<SLP)  and  <  b)  radius  of  30-fcnot  surface  winds 
<R-30>  with  respect  to  average  rate  of 
latitudinal  movement  and  latitude. 
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Figure  32  Climatology  of  the  minimum  surface  pressure 
< SLP>  at  <a>  maximum  intensity  and  at  < b)  radius 
of  30-knot  surface  winds  <R-30>  with  respect  to 
latitude.  Storms  are  at  the  stage  of  maximum 
intensification.  The  solid  line  represents  the 
mean*  and  the  dashed  line  represents  one 
standard  deviation  from  the  mean. 


Ooyaraa's  incompressible  fluid  layer  model  can  be 


rederived  using  isentropic  layers.  This  makes  the 
incorporation  of  thermodynamics  more  straight  forward.  In 
addition*  the  spectral  method  used  to  numerically  solve  the 
governing  equations  allowed  simple  calculation  of  the 
nonlinear  terms. 

The  effects  of  sea  surface  temperature  and  mid-level 
moisture  on  the  development  of  tropical  cyclones 
characterized  by  this  study  are  similar  to  those  predicted 


in  previous 

numerical 

model ing 

research 

( Ooyama* 

1969; 

Sundquist* 

1970; 

Rosenthal  * 

1970* 

1971), 

and 

cl imatological ly  < Gray* 

1968*  1979; 

MacBr ide 

,  1979). 

In 

both  cases*  the  influences  are  monotonic:  increasing  storm 
strength*  size*  and  rate  of  development  with  increasing  sea 
surface  temperature  and/or  mid-level  moisture.  However* 
6torms  are  the  same  size  at  maximum  intensity.  The  greater 
the  sea  surface  temperature  or  mid-level  moisture*  the 
greater  the  amount  of  moisture  available  for  condensation 
and  latent  heating  per  unit  volume  of  ascending  air.  By 
the  CISK  theory  described  by  Ooyama  <1982>*  this  increased 
latent  heat  release  would  produce  a  feedback  which  causes 
faster  and  greater  tropical  cyclone  development. 

Because  the  maximum  heating  is  located  slightly  inside 
of  the  radius  of  maximum  winds*  the  radius  of  maximum  winds 
decreases  with  time  as  the  storm  intensifies  (figs.  7  &  8). 


80 


Because  there  are  greater  amounts  of  heating  occurring 
closer  to  the  center  of  the  storm,  storms  with  moister  mid- 
levels  or  over  warmer  oceans  intensify  faster  and  stronger 
than  storms  with  drier  mid-levels  or  cooler  ocean 
temperatures  (fig.  18).  In  addition,  the  more  intense 
storms  tend  to  deintensify  more  rapidly  because  the  high 
rates  of  heating  required  to  maintain  such  strong  storms  is 
not  possible  with  the  decreasing  V  values  and  because 
surface  drag  is  nonlinear.  Because  heating  of  the  mid¬ 
levels  decreases  the  conditional  instability  of  the 
temperature  lapse  rate,  storms  with  rapid  heating  quickly 
lose  their  ability  to  maintain  cumulus  convection. 

The  effect  of  latitude  is  not  monotonic  as  described 
in  previous  literature  (Yamasaki.  1968).  Yamasaki  only 
studied  three  latitudes  (0°,  4°  and  20°N>  to  describe  the 
effect  of  latitude.  In  this  study,  storms  at  5°N  intensify 
slowly  and  become  very  weak  storms.  Interestingly,  most 
numerical  model  studies  use  20°N  as  the  typical  latitude 
for  tropical  cyclones.  which  is  within  the  bands  of 
cl imatological ly  intense  storms  and  of  numerically 
simulated  strong  tropical  cyclones. 

When  the  latitudinal  variation  of  sea  surface 
temperature  is  neglected,  storms  within  20°  -  30°N  latitude 
intensify  the  most.  yet  do  not  intensify  most  rapidly. 
Storms  at  lower  latitudes  (near  10°N>  arc  very  small  and 
relatively  weak,  and  tend  to  intensify  very  quickly. Whereas 
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when  sea  surface  temperature  or  mid-level  moisture  vary# 
storms  that  intensified  most  also  intensified  the  fastest. 
Latitude  does  not  directly  influence  the  availability  of 
moisture  as  sea  surface  temperature  and  mid-level  moisture 
do,  but  latitude  does  influence  the  positioning  of  the 
maximum  heating  <fig.  33)  and  the  relationship  between  u 
and  v  <  fig.  34).  Sea  surface  temperature  and  mid-level 
moisture  did  not  affect  the  relationship  between  u  and  v; 
the  trend  is  shown  in  Figure  10c. 

Schubert  and  Hack  < 1982)  described  the  effects  of 
inertial  stability  on  the  efficiency  of  heating  within  a 
tropical  cyclone.  The  greater  the  inertial  stability  for  a 
given  heating  field,  the  greater  the  efficiency  of  heating, 
and,  therefore,  the  greater  the  rate  of  intensification. 
Figure  35  shows  that  storms  at  low  latitudes  have  greater 
inertial  stability  at  the  radius  of  maximum  heating  of  the 
storms  in  addition  to  greater  rates  of  heating,  so  these 
storms  intensify  the  fastest.  But  because  the  storm  is 
intensifying  so  guickly,  the  r\  profile  is  being  modified 
more  rapidly  (fig.  36)  and  the  environment  in  the  region  of 
maximum  updraft  and  heating  becomes  less  and  less 
conditionally  unstable,  causing  the  storm  to  stop 


intensifying. 


In  the  case  of  10  N,  the  vertical 


instability,  although  small,  is  large  enough  to  allow  the 
storm  to  maintain  the  maximum  intensity  for  several  days 
(fig.  16).  In  addition,  inertial  stability  is  greatest  at 
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lower  latitudes  so  that  for  even  small  rates  of  heating# 
the  greater  efficiency  of  heating  allows  the  storms  to 
maintain  its  maximum  intensity  which  is  small  compared  with 
storms  at  higher  latitudes.  Also*  surface  drag  is  not  as 
great. 

Tropical  cyclones  at  higher  latitudes  tend  to 
intensify  at  slower  rates  but  intensify  more  than  lower- 
latitude  storms.  A  factor  that  will  cause  northerly  storms 
to  spin  up  more  slowly  is  that  the  radius  of  maximum 
updraft  and  heating  is  farther  from  the  center  of  the  storm 
and  weaker  compared  to  storms  at  lower  latitudes.  The 
weaker  rate  of  heating*  results  in  the  slower  rate  of  spin 
up.  By  the  conservation  of  angular  momentum*  if  the 
forcing  of  the  storm  is  farther  from  the  center  the 
response  will  be  to  produce  a  larger  storm  with  slower 
tangential  winds.  Because  the  heating  is  occurring  further 
from  the  center*  the  conditional  instability  takes  longer 
to  modify.  So  the  northerly  storms  intensify  longer  but  at 
slower  rates.  This  allows  the  poleward  storms  to  reach 
relatively  strong  maximum  intensities. 

When  both  latitudinal  variation  of  sea  surface 
temperature  and  latitude  are  considered*  the  results  of  the 
numerical  simulations  are  very  similar  to  the  climatology 
results  <fig.  21).  Because  storms  at  lower  latitudes  are 
over  waters  with  relatively  warm  temperatures#  the  absence 
of  large*  strong  storms  is  more  a  result  of  latitude.  The 


latitudinal  band  <15°  _  25°N>  where  relatively  strong  larga 

atorma  occur,  is  a  result  of  latitude  affect  and  sea 

surface  temperatures.  North  of  this  band,  the  absence  of 

strong  storms  is  primarily  a  result  of  the  sea  surface 

temperature,  yet  the  larger  storm  size  <i.e.,  radius  of  30- 

knot  winds)  is  a  result  of  latitude.  In  summary,  storm 

size  is  affected  predominately  by  latitude  at  all 

latitudes,  and  the  maximum  storm  intensity  is  affected  by 

latitude  at  lower  latitudes,  latitude  and  88T  at  15°- 
o 

25  N,  and  sea  surface  temperature  at  higher  latitudes. 

The  effect  of  storms  moving  north  from  low  latitudes 
does  not  seem  to  be  cl imatologlcal ly  related  to  how  fast  or 
how  much  a  tropical  storm  develops.  Yet  from  the  numerical 
simulations,  storms  that  are  spinning  up  rapidly  in  lower 
latitudes  would  tend  to  intensify  more  as  the  storms  move 
into  higher  latitudes.  Thus,  as  the  storm  moves  northward, 
the  radius  of  maximum  heating  expands  which  allows  the 
storm  to  intensify  longer,  and  so  stronger.  Further 
studies  of  this  effect  should  be  made  with  a  three¬ 


dimensional  model 
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